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ABSTRACT 

We present 10" to 18" images of four massive clusters of galaxies through the Sunyaev-Zel'dovich 
Effect (SZE). These measurements, made at 90 GHz with the MUSTANG receiver on the Green 
Bank Telescope (GBT), reveal pressure sub-structure to the intra-cluster medium (ICM) in three 
of the four systems. We identify the likely presence of a previously unknown weak shock- front in 
MACS0744+3927. By fitting the Rankine-Hugoniot density jump conditions in a complementary 
SZE/X-ray analysis, we infer a Mach number of A4 = 1.2^0*2 an d a shock- velocity of 1827^195 km s _1 . 
In RXJ1347-1145, we present a new reduction of previously reported data and confirm the presence 
of a south-east SZE enhancement with a significance of 13.9<j when smoothed to 18" resolution. 
This too is likely caused by shock-heated gas produced in a recent merger. In our highest redshift 
system, CL1226+3332, we detect sub-structure at a peak significance of 4.6a in the form of a ridge 
oriented orthogonally to the vector connecting the main mass peak and a sub-clump revealed by weak 
lensing. We also conclude that the gas distribution is elongated in a south-west direction, consistent 
with a previously proposed merger scenario. The SZE image of the cool core cluster Abell 1835 is, 
in contrast, consistent with azimuthally symmetric signal only. This pilot study demonstrates the 
potential of high-resolution SZE images to complement X-ray data and probe the dynamics of galaxy 
clusters. 

Subject headings: galaxies: clusters: individual: RXJ1347. 5-1145, CLJ1226. 9+3332, 
MACS0744. 8+3927, Abell 1835; cosmology: observations; cosmic microwave 
background; GBT 



1. INTRODUCTION 



The Sunyaev-Zel'dovich Effect (SZE) in clusters of 
galaxies arises from inverse Compton scattering of cos- 
mic microwave background (CMB) photons off hot elec- 
trons in the Int ra-Cluster Medium (ICM) (|Sunyaev fc 
Zeldovich||1972[ ). The magnitude of this effect is red- 



shift independent and directly proportional to the line 
of sight integrated pressure of the plasma. At frequen- 
cies < 218 GHz, it is manifested as a decrement in CMB 
intensity. Over the past two decades, measurements of 
the SZE in clusters of galaxies have been used to probe a 
wide range of cosmological and astrophysical questions. 
It has bee n used by dedicated surveys to search for clus- 



ters (e.g., [Hincks et al. 2009; Vanderlind e et al.||2010 



Menanteau et al 



Marriage et al, 



with X-ray data to measure the Hubb 



2010), combinec 
e now (e.g., [Ma- 



|son et al.|200T| iReese et al.|2002| |Bonamente et al.| 2006 ) 

and to der ive physical cluster properties from radial pro' 



files (e.g., LaRoque et al. 2003 Bonamente et al. 2006 
Mroczkowski et al. 20 09]). For revie ws o f the SZE and its 



applications, see Birkinshaw (1999) and Carlstrom et al. 
(2002). 
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Measurements of the SZE at high angular resolution 
are difficult because of the large apertures required. 
Nearly all measurements currently in the literature have 
effective angular resolution larger than ~ V . These an- 
gular scales (corresponding to ~ 365 kpc at z = 0.5) 
are extremely useful for measu ring the bulk sign al out to 
large cluster-centric radii (e.g., Nord et al.|2009[ ) but are 
unable to resolve the smaller scale physical processes in 
the cluster cores. 

High-resolution X-ray imaging from Chandra and 
XMM-Newton in the last decade opened a new window 
to cluster physics. Objects once thought to be spheri- 
cally symmetric and relaxed have been shown to display 
evidence of interesting phenomena which provide insight 
to the complicated dynamics at play in these structures. 
Among these are s hocks and cold-fronts induced by re- 
cent mergers (e.g., Markevitch & Vikhlinin 2007), cavi- 
ties and heating cau sed by AGM interactions (e.g., Mc- 
|Namara et al.||2005| ) and shar p surface brightnes s edges 
caused by gas sloshing (e.g., |ZuHone et al.||2009[ ). 

High-resolution images of the SZE in clusters provide 
a new tool which, when combined with X-ray measure- 
ments, can constrain complicated physics in galaxy clus- 
ters. This is particularly true of the high-redshift Uni- 
verse as the X-ray surface brightness data (proportional 
to the product of the density squared and square root 
of temperature integrated along the line of sight) suffer 
from cosmological dimming. 

The potential of resolv ed SZE was first demonstrated 
by |Komatsu et aL (2001), who used the Nobeyama 45m 
to image RXJ 1347-1145, a massive X-ray luminous clus- 
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ter prev iously thought to be r elaxed and spherically sym- 
metric ( Schindler et al.p997 ). The asymmetry revealed 
by their work was the first indication that the system 
was disturbed, and it is now believed to h ave under- 
gone a recent merger ( Kitayama et al.||2QQ4[ ). Thi s has 
since been con firmed m the SZE by MUSTANG ( |Ma-| 



son et al. 



2010). More recently, other groups h ave m ade 



high resolution SZE images in CL J0152-1347 (Massardi 
et al. 201 0} at ~ 35" resoluti on) and the Bullet Cluster 
( |Savyasachi Malu et al.|2010[ at ~ 30" resolution). 

In this work, we present measurements taken with 
the Multiplexed SQUID/TES Array at Ninety Giga- 
hertz (MUSTANG) receiver on the 100m Robert C. 
Byrd Green Bank Telescope (GBT). The large collect- 
ing area of the GBT combined with the focal plane array 
of bolometers make this system ideal for probing sub- 
structure in clusters through the SZE. 

The clusters MACS0744, RXJ1347, and CL1226 were 
selected because they displayed signs of merger activity 
in previous measurements and therefore were likely to 
contain small-scale features created by merger processes. 
In contrast, Abell 1835 was specifically targeted because 
it was expected to be representative of a typical relaxed 
cluster. All uncertainties quoted in this paper are 68% 
confidence and we assume a cosmology where Ho = 71 
km s" 1 Mpc -1 , Q A = 0.73 and Q M = 0.27. 

2. INSTRUMENT & OBSERVATIONS 

2.1. MUSTANG 

MUSTANG is a focal plane camera with an 8 x 8 array 
of Transition Edge Sensor (TES) bolometers built for the 
Gregorian focus of the 100 m Robert C. Byrd Green Bank 
Telescope (GBT). It has 18.4 GHz of bandwidth centered 
on 90 GHz. The array has a 0.63/ A pixel spacing which 
yields a well sampled instantaneous field of view (FOV) 
of 42" on the sky. More detailed information on the 
in strument can be found on the MUSTANG websit^Uand 
in picker etaT] fl2008| [2009] ). 

2.2. Observations 

Data presented here were obtained during the win- 
ter/spring of 2009 and 2010. The cluster signal was mod- 
ulated predominantly in a "Lissajous daisy" scan pat- 
tern. This strategy was designed to move the telescope 
with high speed (~ 0.5's -1 ) without drastic accelerations 
which can induce feed arm instabilities and pointing wob- 
ble. Faster scan rates move the sky signal to frequencies 
above the low frequency (1//) noise from the atmosphere 
and internal fluctuations. The GBT bore-sight trajec- 
tory during one of these scans is displayed in Figure [IJ 
This observation pattern is centrally weighted and pro- 
duces maps with radially increasing noise levels. While 
the particular scan shown in Figure [I] contains informa- 
tion in a ~ 6' diameter region, only the central ~ 2' 
are well covered. To improve the sky sampling in clus- 
ter cores, each object was mapped with 5 tiled pointing 
centers: one centered on the X-ray surface brightness 
peak and others offset to the north, south, east and west 
by approximately one instantaneous field of view (FOV) 
(~ 42"). Other scan patterns with more uniform s ky cov- 
erage such as the "billiard ball" scan described in |Dicker| 
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Fig. 1. — GBT bore-sight trajectory during a Lissajous daisy 
scan pattern. This observing strategy was used for the majority of 
the observations presented. The movement of the source across the 
sky during the scan has been subtracted. The shaded box indicates 
the instantaneous field of view of MUSTANG 



TABLE 1 

Observation Summary 



Cluster 




Time 
(h) 


Secondary Calibrator 


Abell 1835 


0.25 


3.5 


1415+1320 


Rx J1347.5-1145 


0.45 


3.3 


1337-1257 


MACS J0744.8+3927 


0.69 


5.8 


0824+3916 


CL J1226.9+3332 


0.89 


4.5 


1159+2914 



et al.|fl2009| ), |Cotton et , al. ([2009]) and |Mason et al.| ( |2010 ) 
were used as well. This alternative strategy has the ad- 
vantage of uniform noise across the map, but at the cost 
of slower telescope velocity with sharper turnarounds. 

At the start of each session, out-of-focus (OOF) holog- 
raphy was carried out using a bright (~ 1 Jy) un resolved 
source. This technique, described in detail in |Nikolic| 



7 http: //www.gb.nrao. edu/mustang/ 



et al. (2007), consists of mapping a compact source with 
the GBT secondary in three positions relative to the pri- 
mary: nominally in focus and 3A on either side. An 
automated real-time analysis uses the measured beam 
patterns to fit for phase errors in the telescope aperture. 
Corrections to the active surface of the GBT primary 
are calculated and applied along with pointing and focus 
offsets. 

After every two Lissajous daisy scans (~ 30 minutes 
on source) the beam profile was measured using a nearby 
bright compact quasar. If significant ellipticity or gain 
decrease is detected in the periodic beam measurements, 
the OOF procedure is repeated. These beam maps are 
used in image reconstruction to track fluctuations in the 
telescope gain, atmosphere and pointing offsets. The 
~ 30 minute calibration timescale was chosen as it is 
characteristic of the thermal time constant of the tele- 
scope. The sources used as secondary calibrators for each 
cluster along with the total on-source integration times 
are presented in Table jl] 

Planets were used for absolute flux calibration and 
were mapped at least once per night. The fluxes of these 
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primary calibrators are taken from We iland et al.| ( [2010 ) . 
Several times a night, off-source scans with the telescope 
at rest and the internal calibration lamp (CAL) firing 
with a 0.5 Hz square wave were taken. These are used 
in analysis to fit for the gain of each pixel. The absolute 
flux of the data is calibrated to an accuracy of 15%. 

3. MUSTANG DATA REDUCTION 

A custom imaging algorithm implemented in IDL is 
used to produce maps from the time ordered bolometer 
data. The data is heavily filtered to remove atmospheric 
signal prior to map making. The process is outlined be- 
low: 

1. Gain inhomogeneities across the detector array are 
flat-fielded using the nearest CAL scan. These data 
are also used to identify and mask unresponsive 
pixels (typically 10-15 out of 64). 

2. A template of the atmospheric signal is estimated 
from low frequency fluctuations that are highly 
correlated across the array. This is constructed 
from an average of the time streams from all ac- 
cepted pixels. The model is then low-pass filtered 
in Fourier space to separate the astronomical signal 
on small spatial scales from the atmospheric tem- 
plate. This filtering requires a characteristic fre- 
quency based on the noise properties of the data. 
The template is then subtracted in the time domain 
from each pixel. The effectiveness of this filter re- 
lies on the assumption that the celestial signal is 
not common mode, which is valid only in the limit 
of compact sources. Bulk signals from clusters are 
not well approximated in this assumption. It is 
therefore essential to simulate and quantify the an- 
gular transfer function of the imaging pipeline. 

3. A low-order polynomial is fit and subtracted from 
each time stream. This further removes the long 
timescale fluctuations in the data. 

4. The data all contain a coherent 1.411 Hz signal. 
This is produced by fluctuations in optical load on 
the detectors caused by the thermal cycle of the 
camera's main cryogenic refrigerator. It is well ap- 
proximated by a sinusoid and is removed at this 
stage. 

5. A per-pixel high pass filter is applied in Fourier 
space. This aggressive technique removes all low 
frequency spatial modes from the data indiscrimi- 
nately. A characteristic frequency is defined at this 
stage as well. This will further affect the angular 
scales present in the reconstructed image. 

6. Individual detector weights are computed based on 
the noise characteristics of each detector after the 
processing described above. This is used to create 
an effective exposure time for each pixel on the sky. 

7. The time stream data are then binned on a 2" x 2" 
grid in Right Ascension and Declination. 

The MUSTANG images presented in this work have 
been optimized for peak signal to noise on the compact 
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Fig. 2. — Parameter space optimization for the CL1226+3332 
data. Color scale conveys the peak signal to noise in the exposure 
corrected map convolved with a 10 x/ Gaussian produced with each 
filter parameter combination. 



features of the clusters. As described above, there are 
two selectable filter parameters used in the map maker, 
one for the common mode template and the other for 
the per-pixel high-pass. These selected frequencies cor- 
respond to spatial scales on the sky through the speed 
at which the signal is modulated by the telescope scan 
(usually ~ 0.5's -1 ). The optimal filter for each object 
depends on the intrinsic structure of the source as well 
as the noise properties of the scans used in each obser- 
vation. To determine the optimal filter for each map, 
parameter space is explored systematically by mapping 
each object with varied degrees of filtering. The peak 
signal-to-noise ratio (SNR) as a function of these two 
parameters is displayed in Figure [2j SNR decreases to- 
wards the top right of this figure because too much signal 
is being filtered out. It also decreases towards the lower 
left because too much atmosphere is allowed in the map. 
A single optimal value for each parameter is assumed for 
the entire data set on each cluster. 

The noise in each map is defined by the standard devia- 
tion of all pixels in an off-source region free from obvious 
signal. From extensive Monte- Carlo simulation, we find 
that this calculation provides a good measurement of the 
noise in the map as a whole, provided that it is scaled 
by the square root of the difference in the map weights 
of the areas in question. 

The necessary filtering steps described above result in 
an attenuation of flux in the recovered map. The mag- 
nitude of this attenuation depends strongly on angular 
scale. Typically, all structure larger than ~ 1.5 x the 
instantaneous FOV is removed entirely. The angular 
transfer function for each object mapped is calculated 
using the specific scans and filter parameters selected 
to produce the map. When quantitatively comparing 
model images to observed data, it is essential to apply 
this transfer function to the model before doing so. A 
more detailed description o f the calculation of t he trans- 
fer function is described in Mason et al. (2010). 
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TABLE 2 
Archival Chandra Data 



Cluster 


Time 


Obslds 




(ksec) 




Abell 1835 


222 


495, 496, 6880 






6881, 7370 


MACS J0744.8+3927 


90 


3197, 3585, 6111 


CL J1226.9+3332 


74 


3810, 5014, 932 



4. Chandra DATA REDUCTION 

Archival Chandra data are reduced using CIAO version 
4.2 and calibration database 4.2.0. Starting with the 
level 1 events file, standard corrections are applied along 
with light curve filtering an d other standard p rocessing 
(for reduction details see, |Reese et al. 2010). Images 
are made in full resolution (07 492 pixels) and exposure 
maps are computed at 1 keV. When merging data from 
separate observations, images and exposure maps from 
each data set are combined and a wavelet based source 
detector is used on the combined image and exposure 
map to find and generate a list of potential point sources. 
The list is examined and adjusted by eye and used for our 
point source mask. A summary of the archival Chandra 
data used in this paper is presented in Table [2] 



5. RESULTS 

5.1. MACS J0744+3927 (z = 0.69) 

This massive high-redshift system, found in the MAs- 
sive Cluster Survey [ MACS) of the all sky ROSAT data 
( |Ebeling et al.|2001a| , has ap peared in several studies us- 
ing X-ray and SZE aata (e.g., |LaRoque et al.|2003| [2006; 
Ebeling et al. 2007 ) . Unlike the other clusters in our sam- 
ple, targeted mult i wavele ngth studies of MACS 0744 are 
scarce in the literature. Kartaltepe et al.| ( |2008| ) include 
this object in a red sequence galaxy distribution study 
of a sub-sample of 12 MACS clusters. They note that 
understanding the assembly dynamics of this system is 
made difficult by its complex morphology, which includes 
some evidence of a dense core in the X-ray images, and 
an elongated doubly peaked distribution of red sequence 
galaxies. 

Indications of a hot compon ent are present in the liter- 
ature as well, particularly in |LaRoque et a l. (2003) who 
imaged this system in SZE on arcmmute scales with the 
BIMA/OVRO telescopes. Assuming a gas mass fraction 
within r 5O c0 of fg = - 081 -o.oii> the y obtained a best fit 
SZE temperature of k B T e = 17.9t3° 4 8 keV - The indica- 
tions of hot gas, with a high central density, in an object 
at high-redshift compelled us to include this cluster in 
our sample. It also has the favorable characteristic of 
having no known radio sources in close proximity on the 
sky. 

5.1.1. MUSTANG Data 

The SZE map produced from 5.8 hours of MUSTANG 
data is shown in Figure [3] It consists of a kidney shaped 
ridge ~ 25" long in the north-south direction. From 
east to west, the structure is roughly the width of our 

8 The radius within which the density is 500 times that of the 
critical density of the Universe. 



beam, and thus is not resolved in this direction. The 
curvature of this feature is well described empirically as 
an 80 degree sector of an ellipse with an axial ratio of 
1.25, with the minor axis and center of the observed SZE 
being 12 degrees south of west on the sky. 

5.1.2. Chandra Data 

The Chandra image is shown beside the MUSTANG 
map in Figure [3j It was produced from nearly 90 ksec of 
combined archival data merged from ObsIDs 3197, 3585 
and 6111 and reduced with the method described in sec- 
tion [4j The core of this cluster displays an asymmetric 
X-ray surface brightness morphology with a sharp dis- 
continuity on the western edge. The concave side of the 
SZE peak identified by MUSTANG is aligned concen- 
trically with the convex edge of the surface brightness 
discontinuity in the X-ray. Such an enhancement in the 
SZE in a location offset from the peak in X-ray surface 
brightness requires a significantly heated plasma. 

5.1.3. X-ray Surface Brightness Shock Modeling 

The combined SZE and X-ray image morphology pre- 
sented in Figure [3] is suggestive of a system dominated by 
a merger driven shock- front. Arriving at this conclusion 
based on the existing relatively low SNR X-ray and SZE 
data alone would be quite tenuous; however, the kid- 
ney shaped ridge seen by MUSTANG combined with the 
sharp edge seen by Chandra is difficult to explain with- 
out invoking a shock-heating mechanism. We proceed 
to model the system in the framework of a shock-front 
through a complementary analysis of X-ray and SZE in 
the approach outlined below: 

• The elliptical geometry and location of the shocked 
gas is approximated from the SZE data. 

• This geometry is used to fit a two dimensional X- 
ray surface brightness profile with a model consist- 
ing of three regions: a cold intact core bordered by 
a cold-front, a shock-heated region bordered by a 
shock-front, and a pre-shock region. These corre- 
spond to I, II and III, respectively, in Figure [4j 

• X-ray spectroscopy is performed in each region to 
obtain the plasma temperature. 

• Three dimensional density and pressure models are 
produced from the surface brightness and spectral 
fits. 

• The pressure model is integrated along the line of 
sight to produce a two dimensional Compton yc 
map. 

• A mock SZE image is constructed at the resolu- 
tion and with the angular extent of the MUSTANG 
map, and the model and data are compared. 

We model the X-ray emissivity as a power law, e oc r _p , 
within each region assuming an ellipsoidal geometry with 
two axes in the plane of the sky and one along the line of 
sight (see Appendix |A| for details). The model has 8 pa- 
rameters in total, two characteristic radii, and a normal- 
ization and power law index in each of three regions. We 
perform a Markov chain Monte Carlo (MCMC) analysis 
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Fig. 3. SZE and X-ray images of MACS0744. Left: MUSTANG+GBT SZE at 13'/5 FWHM effective resolution after smoothing. 
Contours are multiples of 0.5cr starting at 3a. Center: Chandra X-ray surface brightness in the cluster core. The image has been smoothed 
with a 1//5 Gaussian. Right: Composite image of Chandra X-ray and MUSTANG SZE. Blue and Green are identical data on different 
logarithmic color scales. Red shows the MUSTANG SZE data. The kidney shaped ridge revealed by MUSTANG is aligned concentrically 
with a sharp surface brightness discontinuity in the Chandra map. 

radii than R s i. 

With both discontinuity radii in hand, the double dis- 
continuity model chains are run with fixed characteristic 
radii. This is enough of a reduction of parameter space to 
produce converged chains. Best fit and 68% confidence 
level uncertainties are shown in Table [3] In this table, the 
parameter / is defined to be the ratio of the normaliza- 
tion of a given region over the normalization in the cool 
intact core (region I). Because the radial dependence of 
the model follows a power law with an exponent less than 
zero, the amplitudes quoted here are normalized at the 
cold-front radius, R s i = 14V 19, to avoid a singularity at 
the origin. Figure [5] shows the X-ray surface brightness 
profile within the fitting region along with the best fit 
model. 

We also ran MCMC fits modeling a constant X-ray 
background in addition to the shock model. It has no 
statistically significant effect on the shock model results. 
This is not surprising as the X-ray background is over an 
order of magnitude down in surface brightness compared 
to the cluster signal at the outermost radius considered 
in the fit. The X-ray background becomes even less im- 
portant towards the inner radii where the cluster signal 
rises. 

To produce Compton yc maps, the three dimensional 
pressure model was nu merica lly integrated along the line 
of sight using Equation | Al 7| out to an elliptical radius of 
60", where the single power law model becomes a poor 
description of the X-ray data. This map is then used to 
produce a predicted SZE image at 90 GHz. After con- 
volving with the GBT beam, the angular transfer func- 
tion of the analysis pipeline is applied to the model in 
Fourier space and compared to the measured SZE data. 
Since the model is only valid in a specified range of an- 
gles about the center of the ellipse, the remaining sky 
w as assumed to be well de scribed by the double j3 model 
of LaRoque et al.| ([2006) and a single temperature of 
8.0 keV. Three model MUSTANG maps were produced 
using this process and are shown alongside the data in 
Figure [6] The model uncertainty is dominated by the 
errors in spectroscopic ksT e in region II. To account for 
this in data comparison we show three model images cor- 
responding to pressure models produced with the best fit 
and the temperature fits to Chandra data at ±lcr. The 



Fig. 4. — Geometry and regions used for elliptical profiles and 
X-ray spectroscopy on MACS0744 overlaid on the Chandra surface 
Brightness image. Green contours are (—4.5, — 5.5)cr SZE decre- 
ment. The three regions correspond to the cool intact core (I), the 
shock heated gas (II) and pre-shock region (III). One X-ray point 
source has been excised from the pre-shock region. The borders of 
the wedge indicate the azimuthal range used in producing radial 
profiles. 



using Poisson statistics for the X-ray data (for analysis 
and statistics details see, e.g., |Reese et al. 2000, [2002} 
Bonamente et al!]|2006 ). Each chain is run for a million 



iterations. Convergence and mixing are checked by run- 
ning two chains and comparing them agains t one another 



dGelman fc Rubinl|1992| |Verde et al 
of burn-in period does not significant 



2003]). The choice 

of burn-in period does not significantly affect the results 
but for concreteness we report results using a burn in of 
10,000 iterations. The model fit is limited to a wedge 
subtending 80° and extending from 10" to 40" from the 
nominal center. This region corresponds to the region 
of interest sugges ted by the SZE and X-ray data as dis- 
cussed in Section IB. 1.21 

Initial attempts to model all 8 parameters at once were 
unsuccessful due to low SNR in these small regions, with 
the chains showing poor convergence. To limit the num- 
ber of free parameters, we implement chains to deter- 
mine the discontinuity radii, R s \ and R S 2, individually 
and then fix those radii. This entails using a single dis- 
continuity model, which has 5 parameters, rather than 8. 
The inner discontinuity radius, R s i, is determined with 
single discontinuity chains using the entire fitting region. 
The outer discontinuity radius, R S 2, is fit with a single 
discontinuity model limiting the fitting region to larger 
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Fig. 5. — Top Left: Chandra X-ray surface brightness elliptical profile (points) and best fit analytical model in MACS0744 (red line). 
Bottom left: Residual of data and model in top left. Blue lines show the best fit characteristic radii for the cold- front and shock-front, 
14'/ 19 and 1&/23 respectively. Top Right: Intrinsic electron number density model produced from the surface brightness fit on the left. 
Bottom Right: Pressure model produced from the above density model and the temperatures derived from Chandra spectroscopy. The 
shaded regions show the uncertainty based on the spectroscopically measured temperatures. Radii in this figure are elliptical and follow 
the conventions described in Appendix [A] 



data model -la model model +la 

> ) I I 

kT (shock) = 13.8keV kT (shock) = 19.7keV kT (shock) = 29.4keV 



-0.4 -0.35 -0.3 -0.25 -0.20 

mJy brn" 1 

FlG. 6. — MUSTANG SZE data and models in MACS0744. Model images were generated by integrating the best fit and la three 
dimensional pressure models from the X-ray along the line of sight and passed through the relavent angular transfer function. Comparison 
to the MUSTANG data shows the excellent agreement with predicted flux scale. A shock temperature closer to the — la value is favored 
by the measured SZE data. 
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Fig. 7. — Multi- wavelength composite image of MACS0744. Green is HST/ACS data in the F814 W band. Red is Chandra X-ray smoothed 
with a 1//5 Gaussian. Blue color-scale and contours show the strong lensing mass reconstruction of [Richard ( 2 011, in prep.| >. White contours 
are the MUSTANG SZE and are identical to those in Figure [3] They are in units of SNR to account tor uneven exposure across the field 
shown here. Galaxy "Gl" is the BCG of the main cluster. The lensing mass reveals a distinct elongation towards the west. Galaxies "G2" 
and "G3" are bright red ellipticals located in the center of a secondary mass peak with no corresponding baryonic emission seen in X-ray. 
The SZE shows no enhancement at this location either; however, the constraint is weaker as the SZE map has large uncertainty at this 
location due to central weighting of scan strategy. "G4" is another bright cluster member which harbors an X-ray point source. It too is 
coincident with a dark matter peak. The presence of peaks in mass distribution with no corresponding baryons is suggestive of a merger 
scenario in which an infalling sub-cluster has passed through the main core, losing its baryons to ram pressure stripping. It is likely that 
the weak shock identified by MUSTANG was produced by one of these events. 

count for Galactic extinction and assume the solar abun- 
dances of |Asplund et~aT1 ([2009]). T he cross sections of 
| Balucinska- Church &; Mc Cammon (1992) with an up 
dated He cross section (|Yan et al.||l998D are used. 

th 



TABLE 3 

Best Fit Parameters for the Shock Model 
in MACS0744 



Region 


/ 


P 


k B T e 
(keV) 


I 

II 
III 


0.480+°^ 
0-406±S;gS| 


0.913t°j™ 
0.986t° o jf 9 

l-15ll°o:SS 


19.71^ 



iated He cross section ( |Yan et al.| | 

statistic, which is similar to tr 



The 



Cash (1979) 



flux scale in the MUSTANG map is completely consistent 
with the X-ray analysis and is suggestive of a tempera- 
ture closer to the low end of the allowed la parameter 
space. 

5.1.4. Chandra Spectroscopy 

Informed by the shock modeling of the Chandra data, 
regions corresponding to the core, shock heated and pre- 
shock regions are constructed and used for spectral ex- 
traction. These correspond to regions I, II and III in 
Figure [4j Since the calibration varies both in time and 
over the ACIS chips, spectra are extracted and response 
files computed for each of the 3 observations individually. 



XSPEC (Arnaud 


1996 Dorman & Arnaud 200 If is 


used to model the 


LCM with a Mekal spectrum flVlewe 


et al.|1985l 19861 Liedahl et al. 1995; Arnaud & Kothen- 


fiug||1985 [Arnaud & Raymond||1992 ). In this fit we ac- 



"cstat" statistic, wnicn is similar to trie 
statistic, is used when modeling the data to properly 
account for low counts. All three spectra are fit simulta- 
neously to the same plasma model with the abundance 
fixed to be 0.3 solar in all cases. The normalizations are 
allowed to float between data sets. The fit is limited 
to photons within the energy range 0.7-7.0 keV. Best 
fit values for the electron temperature and 68% confi- 
dence ranges are summarized in Table [3] Though the 
uncertainty is large in the photon-starved shock heated 
region, it is clear that there is a significant increase in 
temperature in this region compared to the surrounding 
regions. 

5.1.5. Mach Number 

We calculate the Mach number of the shock-front by 
fitting the Rankine-Hugoniot jump conditions. This 
quantity can be obtained independently by fitting the 
jump in density from X-ray surface brightness or in tem- 
perature as measured by spectroscopy. We use the an- 
alytic expressions from |Finoguenov et al.] ( |2010[ ) for the 
Mach number in these two cases 



M p = 



2 P2 
Pi 



7 + 1 -(7-1) 



"ft 



1/2 



(1) 



8 



and 



Mi 



0.05 |r 







1/2 


8^-7 + 

J- 1 


(8|f-7) 2 + 15 





1/2 



, (2) 



where we assume the adiabatic index for a monatomic 
gas 7=| and pi, p2, T\ and T2 are the density and 
temperature before and after the shock. 

The Mach number can also be calculated from the stag- 
nation condition. This relates the ratio of the pressure 
at the edge of the cold- front, P st , over the pressure just 
ahead the shock-front, Pi, to the Mach number through 
the relationship 



Pst 

Pi 



Ml 



7 + 1 



7- 



2Ml 



(3) 



(2002 



as presented in |Sarazin 

We calculate the Macn number for the potential merger 
in MACS0744 using Equations [TJ [2] and [§j The value ob- 
tained from the density jump conditions was calculated 
from the posterior MCMC used in the fit to the X-ray 
surface brightness. This yielded the value M p = 1.21q!2 
where the errors are 1<j and the full discrete probability 
distribution function is shown in Figure [8j The Mach 
number obtained from the relation imposed by the stag- 
nation condition is A4 s t — l-4^o'2 which is m excel- 
lent agreement with the number provided by fitting the 
density jump. The temperature jump conditions at the 
shock yield a higher value, Mr = 2.1+qJ. While this 
measurement suggests a greater shock velocity, the error 
bars are large and it agrees at the 1.3a level with our esti- 
mate from the density jump condition. The flux scale in 
the MUSTANG image suggests the true temperature is 
towards the low end of the Chandra range as is shown in 
Figure 6 The shock velocity in this cluster is 1827^195 
km s _1 assuming the Mach number obtained from the 
density jump conditions. 

5.1.6. Discussion 

This high-redshift system has proved to be an excel- 
lent example of the power of combining resolved SZE 
and X-ray imaging. The high-resolution SZE measure- 
ments reveal a region which is likely the result of a shock. 
Guided by this data, two sharp discontinuities and a 
spectrum consistent with a substantially hotter plasma 
are detected in the low SNR X-ray data. Deeper Chandra 
observations of this cluster will help confirm the presence 
of a shock and more accurately determine its Mach num- 
ber, which for the density jump fit to the current data is 
mildly consistent with a transonic event (A4 = 1). 

Figure [7] shows a composite image of this sy stem in- 
cluding the st rong le nsing mass distribu tion (Richard 
2011, in prep. [ ) (see also | Jones et al. ( |2010[ )). This reveals 
a highly asymmetric elliptical mass distribution elon- 
gated to the west consistent with the r ed sequence mem- 
ber ga l axy distributi o n pre sented in Kartalte pe et al. 
(2008). |Zitrin et al.| ( |2010| ) have also done a mass re- 
construction and independently obtained a similar mass 
distribution. The HST data shown in green contains mul- 
tiple bright red elliptical galaxies with BCG-like charac- 
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Fig. 8. — Mach number in MACS0744 obtained with 3 methods. 
The black histogram shows the discrete PDF calculated from the 
density jump conditions in the MCMC fit to the surface brightness 
distribution. The red area shows the 68% confidence region from 
the calculation of a temperature jump across the shock-front as 
measured from from Chandra spectroscopy. The blue area is the 
68% confidence region obtained by fitting the stagnation condition. 
The solid red and blue lines are the best fit values obtained from the 
temperature jump and stagnation conditions respectively. While 
the stagnation and density jump conditions yield highly consistent 
results, the result from the temperature jump conditions appear to 
be biased high. This is due to a heavy reliance on the spectroscopy 
in the low SNR region II. 

teristics. Galaxy Gl is coincident with the X-ray sur- 
face brightness peak and is assumed to be the BCG of 
the main cluster. Roughly one arcminute to the west of 
the X-ray center, the lensing mass reveals a second peak 
containing the bright red galaxies G2 and G3. While the 
baryon distribution is elongated in the direction of this 
potential sub-cluster, there is no X-ray peak associated 
with it. This is likely explained by ram-pressure strip- 
ping during passage of the sub-cluster through the main 
core. Galaxy G4 is another massive cluster member lo- 
cated west of the main peak. This too has a significant 
dark matter halo with no apparent baryonic peak. The 
presence of multiple peaks in dark matter and galaxy 
density with no accompanying baryonic mass is sugges- 
tive of a merger scenario in which a smaller cluster has 
passed through the main core, stripping it of its baryons 
and producing a shock wave in the ICM. The geometry 
of the westerly elongated multiply peaked dark matter 
distribution is qualitatively suggestive of a merger sce- 
nario in which the shock-heated gas identified by MUS- 
TANG could have been produced. However, an accurate 
interpretation of the merger dynamics requires detailed 
modeling through hydrodynamical simulations. 

5.2. RXJ1347-1145 (z = 045) 

The rich cluster RXJ134 7-1145 is an extremely X- 
ray luminou s galaxy cluster ( Schindler et al.||1997 Allen 



et al. |2QQ2[ ) and has been the object of extensive study 
in SZlk, X-ray, lensing, rad i o and optical spectroscopy 



e.g., [S chindle r~et al.||1997| |Poi ntecout eau et al. 
•" 4 - " 1 onm 1 VUlen et al.||2UU2l ICohen & 



1999 

I Komatsu et al.||2UUl| [Allen et al.||2UU2| I Cohen & kneib 
20021 IKitayama et ai.)|2004| IGitti et al.||20071 |Ota et al. 



[200^ ]Bradac et al.||2008| [Miranda et al.||2008j ). It is 

one of only a tew systems which have been studied with 
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the SZE on sub-arc minute scales thus far and is an 
excellent example of the potential for the SZE to re- 
veal the rich sub-structure exhibited in cluste rs. Initial 
measuremen ts made with ROSAT reported by Schindler 
|et al.| ( p"997[ ) deemed RXJ1347-1145 a relaxed system as 
it showed a round, singly peaked surface brightness mor- 
phology and a strong cool core. SZE observations made 
with the NOBA bolometer camera on the Nobeyama 
45m at 150 GHz (iKomatsu et al. 1120011 iKitayama et al. 



2004) revealed an enhancement to the SZE 20" (170 kpc) 
to the south-east of the cluster center. This asy mme- 
try is now supported by X-ra y and radio data ( |Allen| 



et al. 
TToT{£ 



2002 
bI 



Gitti et al. 2007) and is interpreted as a 
> 20keV) fe ature caused by shock heating in 



a recent merger event ( jKitayama et al. 2004[). This fea- 
ture w as confirmed recently by MUSTANG ( [Mason et al. 



2010). 



The data we present in this p aper are identical to those 
described in |Mason et al. (2010) but are processed with 
the additional per-pixel high-pass Fourier filter described 
in Section [3j By deliberately isolating the high frequency 
spatial modes we are able to increase the signal to noise 
on the small-scale features shown in Figure [9] This comes 
at the expense of removing extended SZE signal. While 
the bulk emission is useful for measuring cluster parame- 
ters like gas fraction and mass, MUSTANG's niche is con- 
straining core sub- structure. S imilar techniques, such as 



un-sharp masking (|Russell et al 



are used in the X- 



the subtraction of a model ( McNamara et al.||2005|) or 

ray to remove bulk emission and examine sub-structure. 

The south-east enhancement is likely due to shock- 
heated gas caused by a merger, but the geometry and 
direction of propagation are not obvious as in the case of 
MACS0744. This makes it difficult to fit the Rankine- 
Hugoniot jump conditions across a discontinuit y in den- 
sity inferred from X-ray surface brightness. Ko matsu] 



et al. (2001) first reported the south-east enhancement 
at a peak decrement of 4.2<j. Figure [9] contains a signal 
to noise map of the MUSTANG data set smoothed to 
comparable resolution to that of the original Nobeyama 
map. This result is now confirmed at a 13.9<r significance 
level in SZE at 18" resolution. 

We find good qualitative agreement of our map to a 
model cluster selected from a suit e of hydrodynamical 
simulations by |ZuHone et al. ( |2009[ ) . This particular sim- 
ulated cluster had recently undergone an off-axis merger 
with a high mass ratio (~ 10 : 1). At the epoch of ob- 
servation, the sub-cluster is moving through the atmo- 
sphere of the main cluster on its first pass of the merger. 
This model also reproduces several other phenomenolog- 
ical features such as a sharp edge in X-ray surface bright- 
ness to the east of the core caused by sloshing of the cold 
gas, and the location of a second bright elliptical galaxy, 
thought to be the BCG of the infalling sub-cluster. 

5.3. CL1 226. 9+3332 (z = 0.89) 

With a mass of (1.4 ± 0.2) xlO 15 M (|Jee & Tyson 



2009) within r2oo, this system is among the largest known 
in the high redshift Univers e. Early measuremen ts of 
the baryons were reported by Ebeling et al. (2001b) who 
identified it in the ROSAT WARPS survey With the 
limited resolution of ROSAT, the cluster was deemed to 
display relaxed morphology. 
Because of its high-redshift, X-ray spectroscopy on this 
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Fig. 9. Top: MUSTANG map of RXJ1347-1145 at 10" resolu- 
tion. Contours are in units of la starting at 3a. Bottom: Signal 
to noise map convolved to 18" resolution. The peak in t he south- 
east quadra nt, originally identified by Nobeyama at 4.2a ( |Komatsu| 
|et al.|2 001), is 13.9<r. 

object is difficult. Initial tempera ture measurements by 
Chandra (Bo namente et al!| |2Q06) indic ated a hot ICM 
(~ 14 keV). This was consistent with iMaughan et al 



(2004) who made previous measurements with XMM 



(~ 12 ke \Q. A more deta i led a nalysis of the ICM prop- 
erties by Maughan et al. ( 2007| ) combined Chandra and 
XMM spectroscopy. They confirmed the hot ICM and 
found an asymmetry in the temperature map with the 
cluster emission south-west of the cluster center hotter 
than ambient. This object has also been mapped in 
the SZE on arcminute scales by the SZA and a strong 
central decrement was measured (Mu chovej et al.| 2007; 
iMroczkowski et al.l[2009l ). 

|Jee & Tyson| ( |2009p mapped the dark matter distri- 
bution ol this system through a weak lensing analysis. 
They found that on large scales the cluster was consis- 
tent with a relaxed morphology but the core was resolved 
into two distinct peaks: the dominant one in close prox- 
imity to the BCG and another ~ 40" to the south west. 
While this sub-clump shows no surface brightness peak 
in either Chandra or XMM data, the location is con- 
sistent with the temperature enhancement reported by 



Maughan et al. (2007) and a secondary peak in the mem- 
ber galaxy density. One possible explanation of this is 
a merger scenario in which a smaller cluster has passed 
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through the dominant core on a south west trajectory 
stripping its baryons and causing shock heating. 

The MUSTANG map is shown in Figure [To] It reveals 
an asymmetric, multiply peaked pressure morphology in 
this high- z system. The most pronounced feature is a 
narrow ridge ~ 20" long located ~ 10" south-west of the 
X-ray peak. A second peak is found in good proxim- 
ity to the X-ray emission which is also coincident with 
the BCG. Also shown in this figure is the X-ray derived 
temperature and pseudopressure (defined as the product 
of the temperature m ap and the square ro ot of surface 
brightness) maps from |Maughan et al.] ( |2007[ ). The Chan- 
dra surface brightness image was produced with 74 ksec 
of archival data taken in Obslds 3180, 5014 and 932. 
There is good qualitative agreement between the two 
data sets, although the small scale features seen by MUS- 
TANG are absent in the Maughan map. This is not un- 
expected as the X-ray derived pressure relies heavily on 
the temperature map which was produced with a variable 
sized aperture. Therefore, adjacent pixels are not inde- 
pendent. This correlation makes the map less sensitive 
to small-scale features. 

Figure [ill shows radial profiles of the X-ray surface 
bright nessToZE and lensing mass distribution. The pro- 
files are centered on the X-ray peak which is coincident 
with the BCG. Each plot shows two profiles, one taken 
from the south-eastern quadrant (red) and the other in 
the south-western quadrant (black). It is clear that all 
data sets are consistent with an asymmetry elongated 
towards the south-west as proposed by the merger sce- 
nario. 

The core of this cluster is compact on the sky due to 
its high redshift. For this reason, the bulk emission is 
likely to contribute non-negligible amounts of flux to the 
MUSTANG map. To quantify the significance of the 
sub-structure, we compare our map t o the best fit spher- 



ically symmetric Nagai et al. (20 0T| model of the SZA 
data as presented by Mroczkowski et al. (2009). Figure 



[I2| shows our map. We assume the spherically symmet- 
ric bulk model is centered on the X-ray peak and take 
the difference between the two maps. The residual map 
figure, contains the peak of the ridge at a 4.6<r level. 

There is a positive unresolved feature in the residual 
map in Figure [12| located 8" northwest of the X-ray peak. 
This unexplained feature could have several interpreta- 
tions. It could simply be a noise artifact. However, it is 
also possible that it is a faint unresolved source. Becaus e 
it is not detected at 30 GHz (Mroczkowski et al.||2009|), 



such a source would require a rising spectrum in the mil- 
limeter as would be expected from a high-redshift, dusty 
star- formi ng galaxy. This g alaxy could be lensed as spec- 
ulated by |Blain et al.] fl2002p and |Lima et al.| ( |2010D and 
similar to the one found in the Bullet cluster by |Wil-[ 



son et al. (2008) and Rex et al 



(20091). Disentangling 



speculation such as this requires the addition of resolved 
millimeter or submillimeter follow-up with different in- 
struments. 

A multiwavelength composite image of this system is 
presented in Figure [13] which includes the weak le nsing 
mass distribution presented in |Jee fc Tyson (2009). The 
northern end of the dominant ridge in the MUSTANG 
image, labelled "B" in this figure, is roughly coincident 
with the lensing mass peak. The orientation of the ridge 
is approximately orthogonal to a vector connecting the 



BCG and secondary lensing peak which is most likely 
the trajectory of the sub-cluster. We posit that the ridge 
is produced by a reservoir of shock-heated gas created 
in the core passage of the sub-cluster, reminisc ent of the 



easte rn pea k in the famous "Bullet Cluster" ( |Markevitch 
|et al.|2002[ ). 

5.4. Abell 1835 (z = 0.25) 

The massive cool core cluster Abell 1835 has proved to 
be an excellent laboratory for studying a range of cluster 
physics. It has bee n used to map the large s cale dark 
matter distribution ( |Clowe fc Sch neider 2002), l ook fo r 
effects of small-scale turbulence (Sanders et al.|[2bl0[), 
search for lensed background submm galaxies (l^vison 
[eTld1|20 00), map the extended radio emission d Govoni 
let al.||2009j and study the central coo l core (e.g., |Pete? 
| son et al.||200l| |Schmidt et al.|[2001[ ). This clust er has 
also been the subject ol extensive SlZE modeling (Reese 



2008 



et aLl[2002l IBenson et al.||2004l |Bonamente et al.| [2006, 
W fHorner et al.||2010p 



Xside from the the central ~ 10" region wh ich displays 
a cavity sys tem excavated by a central AGN (|McNamara| 



et al. 2006), the X-ray morphology is well described by 
a spherically symmetric geometry with no obvious sub- 
structure. This distinguishes it from the rest of our sam- 
ple. The absolute pressur e is extremely high in the core, 
as was demonstrated by Reese et al.| (|2002|) who mea- 



sured a central decrement of — 2.502_ 



.175 



mK at 30GHz. 



However, the MUSTANG map shown in Figure [14] con- 
tains a low signal to noise detection of the SZE. The map 
in this figure has been smoothed to an effective resolu- 
tion of 18" to increase the signal to noise. This figure 
also shows the Chandra image produced from 222 ksec of 
data, merged from Obslds 495, 496, 6880, 6881 and 7370. 
As described in Section|3j the filtering techniques applied 
are optimized to produce high signal to noise maps on 
small-scale structures. Therefore, the lack of high signif- 
icance SZE in the reconstructed image is indicative of a 
featureless, smooth, broad signal. 

Figure 15 shows several pixel histograms taken from 
different areas in the MUSTANG Abell 1835 map. To 
account for uneven exposure, the maps were multiplied 
by the square root of the weight map. This figure con- 
tains a pixel histogram of the central V diameter in the 
MUSTANG map minus an azimuthally averaged version 
of the same map. This is shown alongside histograms 
from a region off-source assumed to contain negligible 
signal and a Gaussian distribution with its a equal to the 
standard deviation of the pixels in the off-source region. 
The residual of the data and azimuthally averaged data 
shows no significant deviation from the noise map and 
is well described by a Gaussian. From this, we conclude 
that the data are consistent with containing negligible 
signal which deviates from spherical symmetry. 

6. SZE FLUX ESTIMATES 

We provide in this section estimates of the SZE flux 
measured by MUSTANG. We quote the integrated flux 
within two radii, #2cr and #3^, which correspond to the 
radii where the mean significance is >2<j per beam, and 
again at >3<r per beam, binned over many beams. We 
compare our estimates to extrapolatio ns of the SZE flux 
reported by |Bonamente et al.] ( |2008[ ) , which was com- 
puted within r25oo using 100 kpc core-cut /3-model fits 



11 





02.0 12:27:00.0 58.0 



02.0 12:27:00.0 58.0 56.0 26:54.0 



a (J2000) 



a (J2000) 



C 30 1 

20.0 
10.0 
33:33:00.0 

50.0 
5 40.0 
30.0 
20.0 
10.0 
32:00.0 
31:50.0 




Maughan et al . Pressure 



02.0 12:27:00.0 58.0 56.0 26:54.0 



D 30 1 

20.0 
10.0 
13:33:00.0 
50.0 
5 40.0 
30.0 
20.0 
10.0 
32:00.0 
31:50.0 




Chandra . 7keV to 7 . OkeV 



02.0 12:27:00.0 58.0 56.0 26:54.0 



20 

18 

16 

14 

12 

10 

8 

6 



+1.0 



+0.8 



+0.6 



+0.4 



+0.2 



-0.2 



-0.4 



-0.6 



6E-8 



5E-8 



4E-8 



3E-8 



2E-8 



1E-8 



B h w S 

oc(J2000) oc(J2000) H 'a c * 

m n o & 

* M +3 N 

£ O 7 

S & U 

Fig. 10. C11226+3332 X-ray and SZE morpholog y. The contours on al l images in this figure are MUSTANG SZE in units of la 
starting at 3cr. Panel A: Temperature distribution from Maughan et al. (2007). Panel B: MUSTANG+GBT SZE image with ll" effective 
resolution. Panel C: X-ray derived pseudopressure map trom [Maughan et al.| ( |2007| . This was produced by taking the product of the 
temperature map and the square root of the surface brightness" Panel JJ: Uhandra surface brightness in the 0.7 keV to 7.0 keV band 
smoothed with a 1//5 Gaussian. 



to the lower resolution 30 GH z OVRO/BIMA data sets. 
We use the |Itoh et ah] ([1998 ) relativistic corrections to 
the SZE flux frequency relation (described in Appendix 
|B| ) when scaling to 90 GHz, assuming the isoth ermal 
temperatures reported by Bonamente et al.| ( |2008 ). Any 
bias due to non-isothermality as well as any discr epancy 
between our t empe ratures and those reported in Bona- 
jmente et al.| ( |2008[ ) leads to < 2% bias in this rescal- 
ing, which is well within the calibration and compact 
radio source contamination uncertainties in both mea- 
surements. 

For A1835, CL1226, and RXJ1347, we bin the flux per 
pixel within a circular region centered on the peak SNR of 
the map. For MACS0744, we choose an elliptical region 
to capture, approximately, the shape of the prominent 
SZE sub-s tructu re (i.e. the shock-heated region reported 
in Section 5.1.1). Uncertainties in the absolute calibra- 
tion are on the order of < 15%, which we include in our 
estimates. In order to account for any compact radio 
sources in these regions, we extrapolate th e flux mea- 



our measurements to obtain estimates of the underlying 
SZE flux reported in Table [5] In the flux estimates pre- 
sented, we have accounted for attenuation by the filtering 
by dividing by the mean amplitude of the angular trans- 
fer function over the Fourier modes between the beam 
scale and the radius to which we are integrating. This 
value is typically ~ 0.7. 

SZE flux provides an estimate, without relying on X- 
ray data, of the thermal energy in the ICM structure 
we see. Table [5] contains the integrated flux and ther- 
mal energy estimates for these four objects calculated 
using the methods described in Appendix [B] We include 
flux estimates at 30 GHz that we re provided by fits to 
the SZA, OVRO, and BIM A data ( |LaRoque etaL|2006 



(Condon et al. 1998) at 



surements from FIR ST (White et al. 1997} and NVSS 

4 GHz and SZA, 



OVRO, and 



BIMA at re 30 GHz ( [LaRoque et aL]|2QQ6| [Mroczkowski 



|et al.|2009| ) to 90 GHz (see Table [I]) using a power law tit 



where the uncertainties are treated in quadrature. The 
predicted fluxes of the radio sources are given in Table [4j 
Variability in source flux is not accounted for although 
the multi-band measurements were obtained at different 
epochs. We are forced to make the assumption of flux 
stability as there are insufficient available time sampled 
radio data to do otherwise. Since we measure SZE flux as 
a decrement, we add these radio source flux estimates to 



Mroczkowski et al. 2009). It is important to note that 
these flux and energy estimates reported here only rep- 
resent structure on scales < 30", which remain after the 
optimal filtering described in Section [3J Contribution to 
the MUSTANG flux from the extended bulk cluster sig- 
nal will be included in this number as well as the two are 
degenerate. Since residual emission from the bulk has 
not been subtracted, these values should be regarded as 
upper limits of the energy contained in the small-scale 
structure. 

While Abell 1835 and RXJ1347 have comparable inte- 
grated flux on large scales, MUSTANG measures a dra- 
matic difference in the cores. This is caused by the large 
amount of signal on small-scales in the sub-structure 
of RXJ1347 and the smooth featureless distribution of 
Abell 1835. 

7. CONCLUSIONS 
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TABLE 4 
Unresolved Radio Sources 



Cluster Field 


Coordinates (J2000) 


a 


Flux(1.4 GHz) 


Flux(30 GHz^ 


Flux(90 GHzg 




a 5 




NVSS/FIRST(mJy) 


(mjy) 


(mJy) 


Abell 1835 


14^01 m 02 s .l +02°52 / 43 // 


2 


31.25±1.57/39.32±1.56 


2.8±0.3 


1.2±0.2/1.1±0.2 


RX J1347.5-1145 


13^47 m 30 s .7 -ll°45 / 08 // 


6 


45.89±1.46/NA 


10.4±0.3 


6.2±0.3 


RX J1347.5-1145 


13^47 m 30 s .l -ll°45 / 30 // 


2 


17.66±3.16/NA 


<0.3 


<0.07 


CL J1226.9+3332 


12^26 m 58 s .2 +33°32 / 48 // 


6 


3.61±0.22/4.34±0.47 


<0.2 


<0.13 



a Coordinates are from FIRST except in RXJ1347 where they come from NVSS. 
b Measured by OVRO, BIMA and the SZA. 

c Extrapolated assuming a power law spectral energy distribution. 



TABLE 5 

Cluster Flux Estimates From MUSTANG 



Cluster Name 




D A 


#2500 






(Gpc) 


(") 


Abell 1835 


0.25 


0.81 


172±5 


RX J1347.5-1145 


0.45 


1.19 


122±| 


MACS J0744.8+3927 


0.69 


1.47 


59±| 


CL J1226. 9+3332 


0.89 
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Fig. 11. Radial profiles of CL1226+3332 from the SZE (top), 
X-ra y surface brigh tness (middle) and lensing mass distribution 
from |Jee fe Tyson] (2009 ) (bottom). Profiles are centered on the 
X-ray peak and are taken from the south-eastern quadrant (red) 
and south-western quadrant (black) . The SZE map was convolved 
with a 10 ;/ Gaussian before averaging. All data sets are consistent 
with an elongation in the south-west direction as proposed by the 
merger scenario 





x 1CT 4 Jy bm" 1 

Fig. 12. Panel A: MUSTANG S ZE map of th e core in CL1226. 
Panel B: The best fit mo del of a |Nagai et a T l |2007) profile to 
SZA data as is presented in Mroczkowski et al. (20U9). The model 
has been passed through the appropriate transfer function. Panel 
C: The residual of Panel A - Panel B. The white contours are 
(—3cr,— Act) which show the significance of the sub-structure not 
accounted for in the azimuthally symmetric model. The cyan con- 
tours are (+3cr,+4cr). The white x shows the location of the X-ray 
peak as measured by Chandra. 



In this paper we have presented high-resolution images 
of the Sunyaev-Zel'dovich effect in four massive galaxy 
clusters produced from MUSTANG observations. Three 
of the four systems probed here display sub-structure 
in the SZE. In the case of MACS0744, we identify a 
likely shock-front propagating with a Mach number of 
A4 = 1.2^q J. The shock- heated kidney shaped feature is 
located between the system's main mass peak and a sec- 
ond peak which shows no evidence of significant baryonic 
mass. In our highest redshift system, CL1226, we find a 
multiply peaked pressure distribution with an asymmet- 
ric morphology. The location and orientation of a ridge 
found in the SZE, along with a south-westerly elongated 
shape, are q ualitatively supporti ve of the merger scenario 
proposed by | Jee fc Tyson| ( |2QQ9| ). We also present a new 
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12:27:00.0 58.0 26:56.0 
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Fig. 13. — Composite image of CL1226. Red shows the Ch andra surface brightn ess in the 0.7keV to 7.0keV band. Blue color scale and 
cyan contours show the surface mass density distribution of |Jee fc Tysonl ( |2009[ ). Contours are linearly spaced in 15 intervals between 
k = 0.25 and hi = 0.59. Green traces the optical emission as measured by the HST/ACS in the F814W band. White contours show 
the MUSTANG measurement in units of 0.5cr starting at 3cr. Location A demarcates the BCG and is coincident with the X-ray surface 
Brightness peak. Location B shows the Dark Matter peak which is coincident with the northern lobe of the SZ ridge revealed by MUSTANG 
imaging. 
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Fig. 14. — Left: MUSTANG SZE image of Abell 1835 smoothed to 18" resolution. Contours are units of 0.5cr starting at 2.5cr. Note the 
central unresolved radio source. Right: Chandra 0.7 keV to 7.0 keV image of Abell 1835 smoothed with a 1.5" Gaussian. Contours on the 
right are identical to those on the left. 
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Fig. 15. — Pixel histograms showing the azimuthal symmetry 
of the MUSTANG map in Abelll835. To account for uneven noise 
distributions caused by non-uniform exposure, histograms are from 
pixels in the map multiplied by the square root of the weight map. 
Green is the histogram of the central 1' diameter with the signal in 
it. The map used was smoothed to 11" effective resolution. Note 
the positive tail caused by the central point source. While the peak 
in SNR of the SZE is moderate in a single beam, the area used for 
this histogram is larger than 20 beams. The black histogram is 
from an area off source believed to contain no signal. Blue is a 
Gaussian distribution. Red shows a histogram over the same area 
as the green line after the subtraction of an azimuthal average of 
the data. Units of the abscissa are multiples of lcr of the Gaussian 
distribution. To account for the different number of pixels in the 
two regions, all histograms are normalized to peak at 1. From this 
analysis we conclude that the MUSTANG signal is consistent with 
containing no significant structure that does not have azimuthal 
symmetry. 



reductio n of the data from o bservations of RXJ1347 pre- 
sented in |Mason et al.|p010| ). This higher signal to noise 
map confirms the previously reported south-east pres- 
sure enhancement at a 13.9<j confidence level. In Abell 
1835 we report a detection consistent with a spherically 
symmetric pressure distribution and no significant sub- 
structure. 



This pilot study has demonstrated the potential of 
high-resolution SZE to identify sub-structures such as 
weak shocks in galaxy clusters. This is particularly true 
of the high-redshift universe where the X-ray data are 
photon starved. A next generation feedhorn-coupled 
TES bolometer array for the GBT is currently in the 
planning stages. With a much larger FOV (4. / 5) and a 
mapping speed 1000 times that of MUSTANG it will be 
able to image a large number of clusters on angular scales 
from 9" to 9' . Other instruments coming on line in the 
next decade, such as ALMA, the LMT, SCUBA2 and 
CCAT will also have high-resolution SZE capabilities. 
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APPENDIX 

A: SHOCK MODEL 

A .1 : Surface Brightness Profiles 

In this work, we measure the density characteristics of a shock-front and cold front in MACS0744 by analyzing the 
elliptical profiles of X-ray surface brightness I(x 1 y) in some observed photon energy band E\ to E2. (In this paper, 
we consider the surface brightness in the 0.7-7 keV band.) Here, we give the analytic expressions for the X-ray surface 
brightness of elliptical X-ray images with discontinuities. The X-ray surface brightness is given by the line of sight 
integral 

1{x > y) = 4tt(i + Zr )v *)<**. ( A1 ) 

where e is the X-ray emissivity integrated over all directions in the emitted energy band E\{\ + z r ) to £"2(1 + z r ). The 
Cartesian coordinates x, y, and z are aligned as shown in Figure [l6j with z being along the line of sight. The cluster 
redshift is z r . The parameter 77 is 4 if I{x,y) is given in energy units, and n = 3 if I{x,y) is in counts units, which is 
generally the case for X-ray observations. 
We fit the data with an analytic expression for the above integral obtained with the following assumptions: 

1. The X-ray emissivity s(x,y,z) is constant on concentric, aligned, similar ellipsoidal surfaces with the geometry 
and conventions described in figure [16} The three principal axes of this elliptical distribution are a, 6, and c. 



2. Two of the principal axes of the distribution (a and b) lie in the plane of the sky, and the third axis (c) lies along 
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the line of sight. We take the x axis of our coordinate system to be parallel to a, and the y axis to be parallel 
to b. The axis given by a or x is along the direction of propagation of the shock and/or cold front. 

3. Between each of the discontinuities, the emissivity varies as a power-law of the radius, e = e Q r~ p . Here, r is the 
scaled elliptical radius r = [(x/a) 2 + (y/b) 2 + (z/c) 2 ] 1 / 2 and p is the power law index. The emissivity changes 
discontinuously at the shock front and/or cold front. 

4. The shock front and/or cold front has rotational symmetry about an axis in the plane of the sky along its 
direction of propagation (c = b). Although we make this assumption in our analysis of the data on MACS0744, 
none of the expressions given below depend on this assumption, and are correct for any c. 

We treat separately each of the regions bounded by one or two discontinuities. In the case of a shock and cold f ront 



there are three separate regions: the pre-shock gas, the shock-heated gas, and the cold front gas. Since equation (Al) 
is linear in e, we can then sum the surface brightnesses of these regions to give the total surface brightness. 

Since the plane of the sky corresponds to a plane of symmetry at z = in this model, the integral for the surface 
brightness can limited to positive z and doubled, giving 

1 f q2 

I(x,y) = —— — / e(x,y,z)dz. (A2) 

V ' y) 27r(l + z r )^ J qi 

The values of q\ > and <?2 > give the extent of the cluster region along the line of sight. The general form for the 
surface brightness for each of the regions obtained with these assumptions after integration is 

'<*■»> -Wr^^"*"*- (A3) 

where we define A to be the two-dimensional scaled elliptical radius 

A{x,y)=^ + £y , (A4) 

and T is the standard Gamma function. The piecewise function (j) takes a form which depends on the complexity of 
the model for a given region of interest. Between the discontinuities, each emission region can be treated as having 
a single outer edge, a single inner edge, or both an inner and outer edge. For example, in MACS0744, the pre-shock 
region has a single inner edge, the cold core has a single outer edge, and the shock-heated region has both and inner 
and outer edge. The total surface brightness is the sum of these three regions. 
For one outer edge, we assume this e dge i s located at r = 1 in three dimensions and at A = 1 in projection. Then, 



the bounds on the integral in equation ( A2 ) are q\ = and 

a/2 



Q2 

and 6 takes the form 



A>1, y ' 



^ = {l-/^-U), A<1 (A6) 

Here, I x (u,v) is the scaled incomplete beta function I x (u,v) = B x (u,v)/B(u,v), B x (u,v) is the incomplete beta 
function, and B(u,v) = T(u)T(v)/T(u + v) is the beta function. Note that very efficient algorithms for calculating 
B(u,v) and I x (u,v) exist a nd can be found a s intrinsic functions on most computer systems. Alternatively, they are 
given in Numerical Recipes ( |Press et al.|1993] ). 
For a single inner edge located at r = K in three dimension and at A = R in projection, the bounds are 

(R 2 -A 2 f \ A < R (A?) 
0, A>R, v ; 



and q2 = oo and we have 



/a2(p-|,|), a<r 



(A8) 

1, A>R. V } 

It is useful to note tha t our expression for a single outer edge is mathematically identical to the expression derived in 



Vikhlinin et al. (2001) but uses the incomplete beta function (which is more convenient numerically) as opposed to 
the hyper-geometric function . 

Finally, for a region with two edges, we will take their locations to be r = 1 in three dimensions and A = 1 in 
projection for the inner edge, and r = R or A = R for the outer edge, where R > 1. The bounds on the integral 
become 

(l-A-f 2 , A<1 (A9) 
0, A>1, y ' 



qi = c 
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and 



The expression for <fi becomes 



q2 = c{ (# -a^\a<r (A10) 




1 1> 

2' 2j 



A<1 

1 < A < R (All) 
A > R. 



A. 2: Density Profiles 



Once we have obtained the power law index p and the normalization e Q by fitting equation |A3 

n. Th 

n e (r) 



to the data, we can 



(A12) 



reconstruct the intrinsic emissivity distribution. This is related to the density distribution n e (r 

' sjr) 1 1/2 

_A(Te,Z)_ 

where A is the X-ray emissivity function which depends on electron temperature T e and abundance Z. 

If XSPEC^Jis used to determine the temperatures in the emission regions, the same models can easily be used to 
determine the value of A. This has the great advantage that the models, temperature, abundances, and instrument 
responses used for the spectral analysis will be completely consistent with those used to determine n e (r). We assume 
here that the model is a single-temperature MEKAL or APEC model. For this purpose, only the shape of the spectrum 
matters, not its normalization, so the region fit in XSPEC need not be identical to the region fit in the surface brightness 
analysis, as long as the spectral shape is assumed to be the same. If the surface brightness / is analyzed in energy 
units, then the procedure is to determine the X-ray flux F of the spectral region in the same band and with the 
same instrument as used to fit the surface brightness. If the surface brightness was corrected for absorption, then the 
absorbing column should first be set to zero. One also needs to record the normalization of the thermal model, which 
is defined as 

10 -u r 

K = Mi+z r m J n ^ dv > (A13) 

where Da is the angular diameter distance to the cluster, n p is the proton number density, and V is the volume of the 
emitting region. Then, the relevant X-ray emissivity function is 

A = F{l + ^ )2 . (A14) 

Here, n e /n p « 1.21 is the ratio of the electron to proton number densities, and is essentially a constant for typical 
cluster temperatures and abundances. 

If the surface brightness is determined in count units (as is typically the case with X-ray observations), then the 
procedure is to set the observed energy band and instrument in XSPEC to the one used for the surface brightness 
measurements, and then type "show" to determine the model countrate CR. Then, the emissivity function is 

CR(l + z r ) 
10^K(n e /n p ) 

A. 3: Pressure and SZE 

With a three dimensional density model obtained through the above procedure and measurements of T e from X-ray 
spectroscopy, one can produce a three dimensional pressure model which can be used to predict the observed SZ flux. 
From the ideal gas law, the electron pressure is simply 

P e (r) = k B n e (r)T e (r) . (A16) 

By integrating this along the line of sight, one can obtain a two dimensional map of the Compton yc parameter 

yc(x,y)= I ^>Ldz. (A17) 
J rn e cf 

Here, ks is Boltzmann's constant, ot and m e are the Thomson cross section and mass of the electron respectively, 
and c\ is the speed of light. 

Assuming that T e (r) is either a constant or is a power-law function of the radius within each region, the electron 



pressure will vary as a power-law of the elliptical radius, and the same analytic expression (equation A3) can be used 
to determine yc(x^y). One simply makes the substitution 

1 - • (A18) 



4tt(1 + z r Y m e cf 

http : / / heasarc . nasa. gov / xanadu / xspec / 
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From a map of yc, it is straightforward to produce a model SZE image. 




Fig. 16. — Elliptical geometry for a single surface brightness edge used in modeling shock fronts. 



B: THERMAL ENERGY FROM SZE 

The surface brightness of a cluster due to the thermal SZE can be expressed, for dimensionless frequency x v = 
hv/k B T cmB where h is Planck's constant, v is frequency, and T CMB is the primary CMB temperature, as the change 
AiszE relative to the primary CMB surface brightness normalization io? as 



A/sze 
Io 



^2 g(xv,T e )k B n e T e dz 

l>eCf J 

lg(xv,T e )P e dz 

letf J 



z -g(x u ,T e ) y c . 



(Bl) 

(B2) 
(B3) 



The primary CMB surfac e brightness normali zation (in units of flux per solid angle) is Io = 2(k B T CMB ) 3 (hci) 2 = 
2.7033 x 10 8 Jy/Sr [see e.g. |Carlstrom et aL] ( |2QQ2[ )]. The factor g(x u , T e ) encapsulates the SZE flux spectral dependence, 
which is a function of electron temperature when relativistic corrections are taken into consideration. In the classical 
physics limit, 



g(x u ) 



(e*» - l) 2 



(B4) 



We integrate the SZE surface brightness in Eq. |Bl| to relate the SZE flux from a region of the sky to the u nde rlying 
electron pressure in the measured ICM feature. The integrated SZE flux is computed (using Eqns. B2, B3, and A17) 



Isze dVt = g(x u )I / yc dQ 



m e cf 



(B5) 



Since dVt = d£/D\(z r ), where dtt is the area integration element, F$ze physically relates to the thermal energy 
content of the gas within a cylindrical region of a cluster (of volume AHAz). The electron pressure P e relates to the total 
pressure P tot by the the elec tron weighting factor ji e « 1.17 (assuming standard abundances) as P tot = (1 + l//i e )P e . 
In terms of the flux (Eq. B5), the thermal energy content is 



th 



3 (l + l//i e )m e c 2 F SZE ^ 2 4 (^) 
2 



For MUSTANG data at 90 GHz, and an assumed k B T e 

E th = \F SZE \D 2 A (z r ) 



a T I g(x u ) 

10 keV and /i e = 1.17, this is: 
3.9 x 10 55 ergs" 



mJy Mpc 2 



(B6) 



(B7) 
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for DA(z r ) in Mpc. In this work, we use the Itoh et al. (1998) relativistic corrections to Eq. B4 
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